Implicit Time Integration Algorithms For High-Order Methods On Unstructured Tetrahedral Grids With p−Multigrid Strategy

نویسنده

  • C. Lacor
چکیده

Computational fluid dynamics has undergone tremendous development as a discipline for three decades. This has been made possible by progresses in many fronts, including numerical algorithms, grid generation and adaptation, as well as the dramatic increase in computer CPU speeds. Moreover, in recent years there is a growing interest in multidisciplinary CFD applications such as aeroacoustics, turbulent combustion and biomedics. These applications are characterized by large scale disparities in the flow fields features, where an adequate and accurate discrete representation is needed. In addition, since CFD is more and more used as an industrial design and analysis tool the applications are often of industrial relevance, requiring unstructured grids for efficient meshing. High accuracy must therefore be achieved on such unstructured grids. Discontinuous Galerkin schemes, Spectral Volume and Spectral Difference schemes are especially suited for these purposes. While relatively good computational efficiency has been attained for the Euler equations, there are still significant challenges remaining for the NavierStokes equations. A major obstacle in achieving such a goal is the geometrical stiffness of the discrete NavierStokes equations caused by the requirement to adequately resolve viscous boundary layers and unsteady propagating vorticies with an economical distribution of grid points. In fact, when spatial numerical methods are combined with classical time integration algorithms, such as explicit Runge-Kutta solvers, there is a restrictive CFL condition, in particular when the cell aspect ratio is high [1]. Moreover, high-order methods are restricted to a small CFL number which can take an excessive amount of CPU-time to reach a state-steady solution with explicit solvers. Hence, with high-order methods, it is possible to achieve low error levels more efficiently than with traditional schemes, but efficient implicit solution approaches are necessary to fully fulfill this potential. In the present contribution an implicit LU-SGS approach is combined for the first time with the Spectral Volume approach and with a p−Multigrid strategy [2]. The spectral volume (SV) method consists of two basic components [3]. One is the data reconstruction, and the other is the (approximate) Riemann solver. In this method in order to perform a high-order polynomial reconstruction, instead of using a (large) stencil of neighboring cells, an unstructured grid cell, (a) Mach contour (b) Residual history Figure 1: Mach contour (a) and convergence history with LU-SGS with Backward Euler difference and optimized explicit R-K smoothers for the inviscid-sine-bump flow (b); 3rd−order SV method; full p−multigrid calculation. called a spectral volume (SV), is partitioned into a structured set of sub-cells, where the cell-averaged solutions are the degrees-of-freedom (DOFs). These DOFs are used to reconstruct a high-order polynomial inside the SV and the numerical flux can be computed with an approximate Riemann solver. The general LU-SGS formulation with Backward Euler difference applied to the spectral volume for the conserved variables Q is indicated as follows ( I ∆t − ∂Rsv ∂Qsv ) ̃ ∆Q sv = Rsv(Q)− I ∆t ∆Q sv (1) where Rsv is the residual. The indices k and k + 1 refer to the previous and current SGS sweeps, respectively, and the symbol (.∗) indicates the most recent solution in the neighbour cells. As can be seen in Table 1 the LU-SGS with Backward Euler difference needs much few iterations and much less CPU-time, being more efficient than the R-K solvers. In the full paper the efficiency of the method will also be discussed for 2D Navier-Stokes cases, showing its efficiency and stability. Computation CPU-time [s] (single grid) CPU-time [s] (full p−multigrid) Optimized Explicit R−K ≈ 258000 (after 105 iterations) ≈ 1328 LU−SGS with BE ≈ 1801 ≈ 414 Table 1: CPU-times for the residual L2-norm to drop 10 orders of magnitude; 3rd−order SV method.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

A p-multigrid spectral difference method with explicit and implicit smoothers on unstructured triangular grids

0045-7930/$ see front matter Published by Elsevier doi:10.1016/j.compfluid.2008.02.004 * Corresponding author. Present address: Departme nautics, Stanford University, 496 Durand Building, Lam United States. Tel.: +1 650 723 6021. E-mail address: [email protected] (C. Liang). The convergence of high-order methods, such as recently developed spectral difference (SD) method, can be accelerated ...

متن کامل

Agglomeration Multigrid and Parallel Processing for the Numerical Solution of the Euler Equations on 2d and 3d Unstructured Grids

Abstract. This paper aims at coupling two known CFD techniques, namely multigrid and parallelization, in an existing Euler equations solver for 2D and 3D unstructured grids. The solver is based on a time-marching formulation for the high-subsonic/transonic flow equations and a pointwise implicit solution algorithm. The gain from the combined use of multigrid and parallelization is the reduction...

متن کامل

Three-Dimensional High-Lift Analysis Using a Parallel Unstructured Multigrid Solver

A directional implicit unstructured agglomeration multigrid solver is ported to shared and distributed memory massively parallel machines using the explicit domain-decomposition and message-passing approach. Because the algorithm operates on local implicit lines in the unstructured mesh, special care is required in partitioning the problem for parallel computing. A weighted partitioning strateg...

متن کامل

An Eecient Cell Vertex Method for Unstructured Tetrahedral Grids

This paper introduces an eecient Euler solver for large tetrahe-dral grids (6:7 10 6 cells). The underlying method is a cell vertex nite volume formulation with a Lax{Wendroo pseudo time{stepping update. This is used as a smoother in an unstructured multigrid method, which requires a sequence of non{nested tetrahedral grids. Fast solutions are achieved by using a a distributed memory parallel c...

متن کامل

A New Implicit Dissipation Term for Solving 3D Euler Equations on Unstructured Grids by GMRES+LU-SGS Scheme

Due to improvements in computational resources, interest has recently increased in using implicit scheme for solving flow equations on 3D unstructured grids. However, most of the implicit schemes produce greater numerical diffusion error than their corresponding explicit schemes. This stems from the fact that in linearizing implicit fluxes, it is conventional to replace the Jacobian matrix in t...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2007